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Abstract 

The throughput achievable in truncated Hybrid ARQ protocol (HARQ) using incremental redun- 
dancy (IR) in analyzed when transmitting over a block-fading channel whose state is unknown at 
the transmitter. We allow the transmission lengths to vary, optimize them efficiently via dynamic 
programming, and show that such a variable-rate HARQ-IR provides gains with respect to a fixed- 
rate transmission in terms of increased throughput and decreased average number of transmissions, 
reducing at the same time the outage probability. 

Index Terms 

Automatic Repeat Request, ARQ, Hybrid ARQ, HARQ, Incremental Redundancy, IR, Block-fading 
channel, Throughput 

I. INTRODUCTION 

Automatic repeat request (ARQ) uses retransmissions to recover data lost due to errors 
inevitable when transmitting over variable and unreliable channels. ARQ is based on the principle 
that the receiver can inform the transmitter about the transmission failure, to which the transmitter 
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responds retransmitting the lost date; ARQ used together with channel coding is known as hybrid 
ARQ (HARQ) HI. HARQ where we limit the number of allowed transmission attempts is known 
as truncated HARQ. 

In this work, we evaluate the throughput achievable in wireless links when using a truncated 
HARQ that conveys incremental, redundancy (IR) in subsequent transmission attempts. For such 
a HARQ-IR system, we use random coding and maximum likelihood decoding assumptions 
of H. We adopt the same simple scenario where each transmission attempt is carried 
out over independently fading channel and we generalize the assumptions of [2] allowing the 
transmission lengths (or - rates) to vary throughout the transmissions attempts. We show how to 
efficiently find the throughput-maximizing rates and we show gains obtained for a finite number 
of transmissions (truncated HARQ). 

The idea of using variable-rate transmissions was already proposed and/or discussed in the 
literature but was not analyzed in the information-theoretic framework of 0, which sets the 
upper bounds on the performance of any practical scheme. For example, a general formulation 
of the problem was provided in which analyzed the infinite number of transmission attempts 
in abstraction of the channel model. The gains of variable-rate transmission over its fixed-rate 
counterpart for the predefined families of code were shown in [JVj [0. In ffTOll ifTTl the 
correlated fading was considered, while lfT2l assumed that the channel stays constant for all 
transmission attempts. The idea of varying the transmission parameters appeared also in lfT3l 
lfl4l lfl5ll Hi, where power was varied on a per transmission-attempt basis. 

We are interested here in the practical case of truncated HARQ when the packet loss (outage) 
cannot be avoided. In such a case the throughput of HARQ may be optimized under constraints 
imposed on the outage probability ifTOl lfT5l or without such constraints lfT3l 0; the latter 
approach is also adopted in this paper. 

In this work we analyze the "conventional" HARQ, i.e., when the return channel can carry 
only one-bit ACK7NACK messages ifTOll . If, on the other hand, we allow the return 
channel to carry more bits, then, the parameters (rate or power) can be adapted using such a 
"rich" or "multi-level" feedback, e.g., lfT6ll Hi- I n me conventional' case, the adaptation is 
not possible but the transmission parameters (rate or power) can be allocated, that is, defined a 
priori for given channel conditions (e.g., the average SNR); this is focus of this work. 

While power adaptation improves the throughput Hi, the power allocation improves the 
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diversity (asymptotic value of the outage for high SNR) but yields significant gains in terms 
of throughput only in the low-SNR range |[T3l . Such conclusions resemble those drawn in the 
context of adaptive modulation and coding IfTTl or in information-theoretic analysis of water- 
filling [fT8l . In this work, interesting in medium-high SNR region, we assume a constant-power 
transmission as the gains obtained when allocating the power are often small |[T3l . 

The objective of this work is thus to evaluate the benefits of constant-power, variable-rate 
transmission for truncated HARQ when compared to the fixed-rate case analyzed in [0 and the 
main contributions are the following: a) we show how to efficiently optimize the rates allocation 
for truncated HARQ with incremental redundancy, and b) we asses the gains of variable-rate 
HARQ over its fixed-rate counterpart, showing that larger throughput, lower outage, and smaller 
average number of transmissions are yield. 

II. System Model 

In the transmission system under study, information bits are separated into packets of equal 
length of N h bits, which are then encoded into codeword of iV s complex symbols x 1 , x 2 , ■ ■ ■ , x^ s 
that are drawn randomly from the zero-mean complex Gaussian distribution with unitary variance. 
The symbols and gathered into K sub-codewords xi,x 2 , . . . , x# whose respective lengths are 
N s ,i, ^V"s,2) • • • j N Bj k. We consider two ways of obtaining the sub-codewords: 

1) A repetition coding (RTC), where the symbols are picked consecutively starting always 
with X\ 

x fc = [xi, . . . ,x Nsk ], N Stk <N s . (1) 

In this way, minfc{iV S) fc} symbols are the same in the transmission attempts 1, . . . , k. 

2) An incremental redundancy (IR) transmission, where each sub-codewords is composed of 
different symbols 

x fc = [x t ' k +i, xt k +N B J with t' k = N s j (2) 

i=i 

This corresponds to puncturing of the codewords x = [x\, . . . , Xn s ] into K distinct sub- 
codewords x fc each of length iV B k = 1, . . . , K, where J2k=i N s ,k — N s and x = [x l5 . . . , x^] 
For convenience, we normalize the values of N Sjk using p k = N Sjk /N h , which has the meaning 
of the redundancy (measured by the number of channel uses per transmitted bit) and satisfy the 

July 3, 2012 DRAFT 



4 



relationship p = N s /N h = J2 k=1 Pk- We define also the rate of each transmission attempt 
R k = 1/pk an d since the rate of the transmission attempts are not the same, we talk about 
variable-rate (VR) transmission, while if p k = px, \/k (or R k = Ri) we obtain the fixed-rate (FR) 
transmission considered before in or Q. 

The ARQ process for each packet starts sending the sub-codeword x x . We assume that 
the feedback (or, return) error-free channel exists, which allows the receiver to send (to the 
transmitter) a one-bit message required by the ARQ process (ACK or NACK). If the packet is 
not decoded correctly^, the NACK message is communicated by the receiver to the transmitter. 
Upon reception of a NACK message, knowing that the first sub-codeword was not decoded 
correctly, the transmitter sends a sub-codeword x 2 composed of iV s 2 symbols. After unsuccessful 
decoding, another NACK message is generated to which the transmitter responds sending the 
codeword x 3 . This continues till the maximum allowed number of transmission attempts K is 
reached (truncated HARQ) or until an ACK message, denoting a successful decoding, is received. 

In a particular case of p k = pi, the sub-codewords have the same length/rate and we recover 
the retransmission schemes analyzed in 0. 

The channel remains constant during transmission of the A;th sub-codeword k — 1, . . . , K and 
the received signal is given by 

Yk = y/^^k + z fc (3) 

where z k is the vector of zero-mean complex, unitary-variance uncorrected Gaussian variables 
(modelling noise). The signal-to-noise ratio (SNR) j k defines the channel state information (CSI) 
which is perfectly known/estimated at the receiver, but unknown to the transmitter. SNR does 
not change during the transmission of the sub-codeword but varies independently from one sub- 
codeword to another. This corresponds to a practical scenario where subsequent sub-codewords 
are not sent in adjacent time instants and, being sufficiently well separated, the realizations of 
the SNR become — to all practical extent — independent. 

'The receiver can determine if the decoding error occurs using an outer error-check code which causes the transmission 
overhead which we neglect for simplicity of the analysis. 
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The channel gains ^7 are Nakagami-m distributed, so the SNR is characterized by the gamma 
function (PDF) 




where 7 is the average SNR. The cumulative density function of SNR is thus given by 



with r(m, 7) = x m ~ 1 e~ x dx and T(m) = T(m, 00) are, respectively, the incomplete 

gamma function and the gamma function. 

The coding scheme is revealed to the transmitter, which in the kth transmission implements 
a maximum likelihood decoding using the observations y fc = [y 1; . . . , y k ]. 

The system-level implementation of the variable-rate HARQ described above deserves some 
comments. Namely, we may assume that each transmission contains only one sub-codeword in 
which case the duration of transmission attempts must vary. This might be a valid approach for a 
single-user communication where the transmitter and the receiver can negotiate the transmission 
time for each sub-codeword. On the other hand, it may be a questionable strategy in multi-user 
communications, where sharing the requirement for a variable-rate transmission with all the 
users is not practical. It might be possible to assign the resources (time) independently of the 
varying transmission length but it would lead to the bandwidth loss (sub-codewords shorter than 
the assigned transmission time slot) or to collisions (sub-codewords longer than the available 
time). 

To avoid such a conceptual difficulty, we assume that the sub-codewords corresponding 
to different packets are gathered in frames that have the duration of Np symbols. Such an 
assumption, also used in |fl9l , ll20l allows us to deal with variable-rate codewords to fill up the 
frame and corresponds to TDMA-type communication, where users are provided with a fixed 
transmission time (frame). This is shown schematically in Fig. [Q We can easily see that the 
relative loss due to variable length of the sub-codewords can be made arbitrarily small, provided 
the number of packets in each frame is sufficiently large. 

III. Achievable throughput 

The definition of the throughput we use here follows 0; according to the reward-renewal 
theorem J2TJ it is the ratio between the expected number of correctly received bits (after up to 
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K transmissions) and the expected number of channel uses iV s required by the HARQ protocol 
to deliver the packet (in up to K transmission attempts). 

We denote by NACK fe , the event of decoding failure in the A;-th transmission and by = 
Pr{NACKi, . . . , NACKjt_i, NACK^} - the probability of decoding failure after k transmission 
attempts. The throughput can be then expressed as flU 

Vk(pi, ■ ■ ■ , Pk) = t^J K f • ( 6 ) 

Pi + l^k=2 Jk-lPk 

which generalizes the results of [|2|j to the case of transmission with variable sub-codewords' 
lengths. Note that fx has the meaning of "HARQ outage", that is, the probability of loosing the 
data packet after the HARQ process is terminated. 

The formulation © is entirely general and depends only on the model of the channel and on the 
coding/decoding scheme. For example, it was used in [|9]| for convolutionally coded transmission 
while J2l used it in independently block-fading channel assuming that capacity-achieving codes 
are available but under constraint p k = p x . Here, we remove this constraint but still follow the 
approach of [2] that has the virtue of providing limits to any practical coding/decoding scheme. 
We thus assume that the coding/decoding scheme is "capacity-achieving" in the sense that the 
transmission is successful if the effective transmission rate is not greater than the accumulated 
mutual information between the sent and the received signals |^ This assumption as well as the 
way the transmitter/receiver deal with the retransmissions will affect the variables fk used in 
©. Namely, three HARQ schemes are considered: 

A. HARQ-I 

In HARQ type-I (HARQ-I), after k transmissions, only the most recent received sub-codeword 
is used for decoding and others are discarded (in [2] this scheme was denoted as ALO). In such a 
case, the decoding failures are independent of each others and the probability of losing a packet 
after k transmissions is calculated as 

k k 

fa = Y[Pi{C( ll )p l < 1} = Y[u( Pl ) (7) 

1=1 1=1 

2 The existence of the codes satisfying this criterion when JVb — > oo is discussed, e.g., in 1221 1161 . 
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where where C(j) = log 2 (l + 7) is the average mutual information (per channel use) when 
transmitting with SNR 7 and u(p) = F(2 1 / p — 1; m) is the probability of outage (after a single 
transmission) when transmitting with redundancy p. 
The throughput of HARQ-I is then given by 

Pl + £fc= 2 P*IL=l "(Pi) 

and the optimal throughput is denoted as r\^ K = max pir .. tPK r}i tK (pi, ■ ■ ■ , Pk)- 

Proposition 1: The maximal throughput of HARQ-I fji t K is independent of K, i.e., fj^K = 
r)i = 771,1 (pi) where p\ = arg p max x ~ v ^ , and is yield with fixed-rate HARQ (FR-HARQ-I) 
Pi,i = pi,l = I, ■ ■ ■ ,K. 



1 ^(p ) whcr^ the ec[u°^ 
in ©, which yields the following inequality 



Proof: Since 171,1 > — y^-, where the equality hold only for p = p\, we can use p k > 



1 - nJLi u (pk) 



Vi,k{pi, ■ ■ ■ , Pk) <Vi,i 



Vi,i 



1 - v{px) + Efc=2(! - U (Pk)) Utl u (Pi) 



1 - n fc =i^(pfc) 



1 - v(pi) + Ef=i 1 nf=i "(pi) - Ef= 2 nti v (pi) 



Vi,k(pi, ■ ■ ■ , Pk) < m,i = Vi (9) 

thus 571 is the maximum throughout of VR-HARQ-I, achievable only if p^ — p\, k — 1, . . . , K. 

According to Proposition [T] the fixed-rate HARQ-I is optimal so the same sub-codeword may 
be used for each transmission and the transmitter can apply the RTC transmission scheme defined 
in Sec. M 

Proposition 1 that is valid for any K may be seen as a generalization of Corrolary 1 in 
valid for K — > 00. 

B. HARQ-IR 

In incremental redundancy HARQ (HARQ-IR) the transmitted sub-codewords are obtained 
according to IR principle described in Sec. fll] and the decoding fails in the k-th transmission 
attempt if the accumulated mutual information is lower than the transmission rate, which yields 
the following condition |fT6ll 

k 

fm, k = Pr{j2C(li)pi<l}. (10) 
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where 72 is the SNR during /th transmission attempt. 

To calculate /ir^ we may proceed as suggested in [2] introducing random variable Vi = C(ji)- 
pi, I — 1, . . . , k whose PDF can be obtained by definition as gi(x) = \n(2)-p(2 x / pl — 1; m)2 x l pl / ' p\. 
This is what will be called the "exact" calculation. 

Alternatively, we may approximate V\ by a Gaussian variable, 0), i.e., 

( \ ~ 1 \ 1 ( ( x ~ c mpi) 2 

gl\X) « gfc(x) = -= 

V 2\ 

where 



i(x) « g fe (x) = -= exp — ^ ) (11) 

V2tv pia m V ^cr^ 



m 



C m = / C , (7)p(7;^)d 7 (12) 

POO 

a 2 m = / C 2 ( 7 )p(7;m)d7-^ (13) 



are, respectively the mean of C(j) (i.e., the ergodic capacity), and the variance of C(j). 

Since v t , I = 1, . . . , k are independent, /i Rjfc = Pr |X)?=i ^ < l| = g fe (x)dx, where g fc (x) 
is a convolution of g^(x). The latter must be calculated numerically, e.g., via direct/inverse 
Fourier transform if the exact form of gi(x) is used, while, applying (fTTI) we obtain a closed- 
form approximation of (flOl) 

/nLfc « /m,fc = Q (t^b 1 ^) , (14) 



A' 



where Q(x) = ^ £° exp(-t 2 /2)dt, £ = £, X k = Ef=i rf, n = VEti P'z , and = p ■ C. 

Proposition 2: Denoting by i]m,K{pi, ■ ■ ■ , Pa:) the approximation of the throughput obtained 
using ([141) in ©, the following inequality holds 

?7ir,a-(pi, • • • , PA") > ?7iR,A-(pi, • • • , Pa-) = C- — g m,i ^ (15) 

Pi + 2^fc=2 Pfc ' /m,fc-i 

where / IR)fc = Q(f (1 - 1/X fe )). 

Proof: The obvious relationship Y" fc 2 < X| used in (fl4)) yields /ir^ < fm,k- From this 
inequality and knowing that the throughput decreases monothonicaly with fk (if pk are kept 
constant), we immediately obtain the lower bound (fl"5l) . 

The bound (fl"5l) will be useful to optimize the throughput in Sec. [IV] 
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C. HARQ-CHASE 

Instead of discarding the packets that were not decoded correctly (as done in HARQ-I), the 
receiver should take advantage of all received packets and if RTC is employed, the received 
signals should be weighted by the corresponding SNR and added up. This is known as maximum 
ratio combining (MRC) or Chase combining ll23l . Then, the decoder used the following signal 

k 

Yk = ^2Vlry , i (16) 

1=1 

where y' k = [y k , 0, 0, . . . , 0] are zero-padded version of the received signal y k . The padding is 
used to make the notation compact and may be seen as an operation carried out at the receiver 
thus it does not affects the throughput. 

Although, in the case of a fixed-rate HARQ, it was shown to introduce little gain over HARQ-I 
ED, lfl3l . Chase-combining is the most the receiver can do when RTC is implemented at the 
transmitter so we deal with this case for completeness of our analysis. 

Calculation of the decoding failure probability is slightly more involved in this case. 

First, for convenience, we reorder the variables 71 , . . . , 7* so that the corresponding sub- 
codewords lengths' after reordering are non-decreasing p Kl < p K2 < . . . < p Kk , where rzx, . . . , K k 
is a permutation of 1, . . . , k. We emphasize that the reordering is merely a concept simplifying 
the analysis and not an a priori constraints on the values p k . 

As illustrated in Fig. [2l thanks to the reordering we are able to identify k "chunks" of the 
sub-codewords 

e M = [xi, . . . ,x NaiKi ), 

£k,l = [ aJ iV s , )tl _ 1 +l: • • • i X N BtKl ], I = 2, . . . , k, 

each with normalized redundancy pi = p Kl — p Kl l (we set p Ko =0), such that all symbols in 
the chunk e k j were transmitted in the transmissions attempts indexed with k,i, ki + i, . . . , n k . 
After simple algebra, the combining of the received signals (fT6l) yields 

fk = [y\i,yfc,2, • • • , y*,fc] (17) 



Yk,i — v7M 



lk,i ■ £k,i + Sbj 



(18) 
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where 

k 

f=i 

is the equivalent SNR for the chunk e k j and £ fc t is a zero-mean, unitary-variance Gaussian vector 
modelling "equivalent noise" affecting the chunk e kji . 

Since the symbols in the chunks e^i are mutually independent, the parts y^i of the received 
signal yk may be seen as the result of transmission of the chunks over the channel with SNR 
jkj. Consequently, Chase combining may be seen as a form of IR transmission with redundancy 

pi and the probability of the decoding failure is given by 

k 

/GH,fc = Pr{X;C f (7 fc ,,)A<l} (20) 
i=i 

that, in the case of a fixed-rate transmission, boils down to the formula shown in 0. Namely, 
since in fixed-rate HARQ-CHASE pi — pi, pi — 0, 1 — 2, . . . , k, and 7^1 = Yli=i 7fc> then the 
decoding failure is calculated in a closed form 

fcu,k = Pr j^E^) • Pi < l| = F(2 1/Pl - 1; m ■ k). (21) 

While d20l) resembles (flOl) . the equivalent SNRs 7^^ appearing in (|2Q|) are not independent 
(unlike in the case of HARQ-IR), so the approach of Sec. IIII-B L based on the convolution of the 
individual PDFs cannot be applied and a multidimensional integration over 71, . . . , 7^ is required 

/cH,fe = / J p(7f, m)djt . . . d7fc (22) 

where V k = {71, . . . , 7* : E/Li C{lk,l)Pi < 1}, so 

fcn,k= / p(7i; m)d7i . . . / p(7 fc _i;m)d7 fc _i / p(7 fc ; m)d7 fc , (23) 
^0 Jo Jo 

where the integration limit for the SNR 7/ depends on the values taken by the SNRs "fi+i, • • • , 7fe 



zi = zj(7 I+1 , ...,lk) 



,^(l-E/=«+i ^[«?(7/)-C(7/+i)]) 



7J+0- (24) 



To implement (1231) we used the Gauss-Laguerre formulae with 10 (for m = 1,2) or 40 (for 
m = |) points in each of fc dimensions of 

The multidimensional calculation was particularly computationally-intensive for K > 4 and 
the results do not seem very relevant beyond this point as virtually all improvement is due to 
the second transmission attempt. 
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D. Limiting cases 

We know from that for a fixed-rate HARQ-IR 

Vm,K > C m (25) 

K— >oo 

where C m is the ergodic capacity, defined in (1251) . 

We also known that, for a given set of pi, . . . , p K , the relationship > /cH,fc > /nt,fc holds 
for all /c 0, and since, for the given pi, . . . , p#, the throughput r\ K © monotonically decreases 
when f k increases, we conclude that for K < oo 

Vi,i = Vi,k < Vck,k < Vir,k < C. (26) 

Thus, the throughput of HARQ schemes with fixed-power transmission, operating without knowl- 
edge of instantaneous SNR, is lower-bounded by ffi defined in Sec. IIII- Al and upper-bounded by 
the ergodic capacity C. 

The limiting case K — > oo is also interesting since, as stated in [[51 Lemma 1], when the 
receiver does not discard packets (as it is the case for HARQ-IR and HARQ-CHASE), the 
optimal redundancy sequence must be non-increasing, i.e., p IR fe > Pir^+i and pcn,k > PcH,fe+i- 

IV. Optimization 

The "design" of the HARQ scheme consists in the maximization of the throughput over 
the redundancy values p 1 , . . . , p K . In the case of FR-HARQ the exhaustive search over one- 
dimensional space is relatively simple. On the other hand, the solutions for VR-HARQ-IR and 
VR-HARQ-CHASE are more difficult to find as their require a multidimensional optimization. 

To maximize © we might use a gradient-based method but the initialization of the variables 
is critical to ensure rapid convergence and to avoid getting trapped far from the global optimum 
(both - not guaranteed in non-concave functions we deal with, cf. [|24l Fig. 1]) , so we used this 
approach only in VR-HARQ-CHASE where various initializations were tested and the solutions 
were compared to the random initializations. This was tedious but feasible as it was done only 
for K < 4. 

In case of VR-HARQ-IR, different approach was adopted: instead of maximizing the through- 
put ij]R t K we maximizing the lower bound (fl"5l) . The problem is greatly simplified since each 



July 3, 2012 



DRAFT 



12 



term / IR)fc depends uniquely on X k = Ya=i p'i an d the optimization may be written as 

max 77i R (pi, . . . , p K ) = max (27) 
pi,-,PK x V K {X) 

where / fe (X) = /m,* = Q(£s/fc(l - 1/X)) and 

V fc (X)= mm pi + V^^.!) (28) 

fe-i 

= min min p' x + V p'Ji^Xi-i) + p' k f k -i(X - p k ) (29) 
o<p' fc <x pi,...y fc _ i: ^ 

Htl p'i= x -p'u 

= min ^(X-^+p/^^X-p). (30) 

0<p<X 

For a given X, the minimization in (l30l) is done over one variable (p = pf.) provided the 
results of the minimization Vk-i(X) are known for all arguments X. That is, first we solve 
V 2 (X) = mm p {X - p + pf x {X - p)}, next V 3 (X) = mm p {V 2 (X - p) + pf x {X - p)}, etc. 
This recursive formulation is characteristic of the so-called dynamic programming (DP) E51 
whose application for throughput optimization was already suggested in 0. The direct benefit 
is that the optimization (1271) over X-dimensions is reduced to K, one-dimensional functional 
optimizations, which greatly simplifies the implementation. 

The function Vk(X) is not obtained in the closed-form, so we discretized X using 50-100 
points over the domain X E (0, k), where the bounding of X by k is not restrictive and comes 
from the heuristic observation that p' k < 1, i.e., each rate = 1/p^ is greater than the ergodic 
capacity C. 

The optimization results are stored as pk(X) = argmin p Vk-i(X — p) + pfk-i(X — p), so once 
the functions Vk(X) are obtained, we can recover the solution p' k that maximizes the bound: 

p' k = p k (X k ) (31) 

where X K = arg x max ±^ and X k _, = X k - p k (X k ). 

We note that while the approximate expressions for /i R fc and ^ir^ are used in DP optimization, 
the throughput values we show in the following are based on the exact calculation of fm tk . We 
also verified that using the DP-based results as the initialization to the gradient-based optimization 
yields practically the same values of the throughput as those we show. 
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V. Numerical Results 

The optimized throughput of fixed- and variable-rate HARQ-IR is shown in Fig. [3] for K = 
2,4,8 for Rayleigh fading channel (i.e. with m = 1), where the gain due to variable-rate 
transmission is particularly notable for HARQ-IR while it is very slight when considering HARQ- 
CHASE, which at best (with K = 4) equals the performance of FR-HARQ-IR with K = 2. 

The gain in terms of throughput offered by VR-HARQ-IR is particularly clear for K = 2 and 
to complement the results of Fig. [3l we evaluate and show in Fig. [4] the "residual throughput" 

X = 1 - B (32) 

C 

i.e., the relative gap between the throughput attained with up to K transmissions and the 
maximum achievable throughput (ergodic capacity). The relative gain of the VR-HARQ with 
respect to FR-HARQ remain roughly constant for all K but of course the absolute difference 
diminishes - as expected - with K since, asymptotically both schemes are equivalent. These 
gains are also more notable when increasing m. The "saturation" of the throughput of HARQ- 
CHASE scheme is also clearly shown. 

In Fig. [5] we show the normalized redundancy p' k = pk ■ C directly proportional to the 
subcodewords' lengths N Bk (inversely proportional to the transmission rates R k ). We observe 
that the first transmission attempt of VR-HARQ-IR is carried out with the rate R\ = 1/pi close 
to C, while the rates of subsequent transmissions increase (i.e., the subcodewords are shorter) 
and decrease again for k approaching K. This relationship holds for all 7 and m and may be 
observed in the IR and CHASE schemes. In Fig. [6l we reproduce similar results for VR-HARQ- 
IR and FR-HARQ-IR but for different values of K. The same "profile" of the redundancy is 
obtained for all K and we may also appreciate that the values of p' k are decreasing with K, which 
is consistent with the optimal behaviour for K — > 00, when the optimal sequence of p k should 
be non-increasing [[51 Lemma 1]. For the FR-HARQ-IR, we observe that p[ decreases with K. 
Recall that, according to the proof in [2, Appendix C], when K — > 00 the throughput-maximizing 
redundancy p' k = p\C should tend to j^. 

The decreasing-increasing behaviour of the values p' k can be interpreted from © combining the 
results of Fig. [6] with those in Fig. [7] showing the values of the decoding failure fk,k = 1, . . . , K. 
Namely, as we strive to make r] K approach closely C, from © we conclude that redundancy/rate 
should be allocated so that C ■ (pi + J2k=2 Pk ' /fc-i) — Pi + Sfc=2 Pk ' fk-i grown to be as close 
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as possible to 1 — fx « 1. Immediately we conclude that we have to use p[ < 1 (transmission 
rate Ri > C) but the behaviour of optimal values p' k ,k>l depends on how the values fk evolve 
with p' k . 

In the particular case of K — > oo, as long as the receiver "accumulates" the redundancy, the 
optimal values p' k should be decreasing with k ll5lR Thus, the fact that p' k increases with k (here: 
for k > 2) is due to the truncation (finite K) and reflects the fact that not only the denominator 
of © should be minimized but also we have to guarantee that the value of fx remains small. 

Also, since fk decreases mach faster in HARQ-IR than it does in HARQ-CHASE (due to lack 
of additional information coded symbols conveyed in the subsequent transmission attempts), the 
optimal values of p' k ,k > 1 can be smaller for VR- HARQ-IR than they are for VR- HARQ- 
CHASE so, as shown in Fig. [5] the variation of the redundancy is less pronounced. 

In Fig. |7] we can also observe that for sufficiently large K (K > 4), the probability of outage 
fm,K m VR-HARQ-IR is smaller than in the case of FR-HARQ-IR. For other values of m and 7 
the same property was consistently observed which is another clear advantage of VR-HARQ-IR 
over FR-HARQ-IR. 

Another consequence of using short sub-codewords for all transmission attempts in FR-HARQ- 
IR is that the mutual information accumulates "slowly" with the retransmissions. Consequently, 
the failures in the initial transmissions occur more likely than in the VR-HARQ-IR, where the 
first transmission is done with the rate Ri close to C. This impacts the average number of 
transmissions which we calculate as 



K-l 

K avg = 1 + J2fk (33) 

k=l 

and show in Fig. [8l 

We can appreciate that when the number of transmission K grows, the average number of 
transmissions i^ avg increases as well but is significantly greater for fixed-rate HARQ-IR: it 
practically doubles for K — 8 and 7 = 30dB. Since the average number of transmissions is 
related to the packet delivery delay (as retransmission can be done only in separate frames), VR- 
HARQ-IR -besides the increased throughput- offers an additional advantage over FR-HARQ-IR. 

3 Remember that for HARQ-I, i.e., when the receiver discards the redundancy of past transmission attempts, the optimal 
solution is pk = pi 
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Figure 1. Example of the structure of three frames sent over channels with corresponding SNRs 71, 72, and 73 when delivering 
data packets denoted by p, I — 1, . . . , 9. The subcodewords having different lengths are identified with different colors and 
patterns. The first frame is filled up with subcodewords of length 7V s ,i (thus, in our example, Nf = 6JV B ,i) corresponding to 
the packets Pi — Pa. When transmitting this frame with SNR 71, we assume C (71)^1 < 1, consequently, the decoder fails to 
decode the message in the packets Pi — Pa and a NACK messages are sent to the transmitter. The next frame contains thus 
six subcodewords of length JV 8j 2 each carrying the redundancy for the undelivered packets and since, here, N B ,i > N s ,2, the 
"empty" space is filled with two subcodewords of the length jV s ,i corresponding of the packets P7 and P$ that are ready for 
transmission. None of the packets is decoded after the transmission of the second frame so, again, six sub-codewords of length 
7V Si 3, corresponding to the packets Pi — P3 are sent as well as the sub-codewords of length N s ,2 corresponding to the packets 
P7 and Pg. The residual time is filled with the sub-codeword corresponding to the packet P9. Note, that the relative loss due 
to unshaded/unfilled space can be made arbitrarily small loading the frame with many sub-codewords. 



VI. Conclusions 

In this paper we have analyzed HARQ with incremental redundancy (HARQ-IR) for transmis- 
sions over block-fading channels. We have proposed an efficient method to allocate the optimal 
ratesand have demonstrated that the variable-rate HARQ-IR provides gains over the fixed-rate 
HARQ-IR in terms of increased throughput, lower outage, and decreased average number of 
transmissions. 
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Figure 2. In RTC, the parts of the transmitted sub-codewords that contain the same symbols are identified as "chunks" ek,i', 
symbols in each chunk experience the same equivalent SNR 7^; = Ylf=i 7«/- 
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Figure 5. Throughput-maximizing normalized redundancy p' k = p • C, k = 1, . . . , K (K = 4) for VR-HARQ-IR (solid, green 
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Figure 7. Outage values k = 1, . . . , K for VR-HARQ-IR (solid, green line) and FR-HARQ-IR (dashed, red line); m = 1, 
7=10dB. 
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